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We use the coupled cluster method to investigate the ground-state (GS) properties of the frus- 
trated spin-^ J1-J2-J3 model on the honeycomb lattice, with nearest-neighbor exchange coupling 
Ji plus next-nearest-neighbor (J2) and next-next-nearest-neighbor ( J3) exchanges of equal strength. 
In particular we find a direct first-order phase transition between the Neel-ordered antiferromagnetic 
phase and the ferromagnetic phase at a value J2/J1 = —1.17 ± 0.01 when Ji > 0, compared to the 
corresponding classical value of —1. We find no evidence for any intermediate phase. From this and 
our previous COM studies of the model we present its full zero-temperature GS phase diagram. 
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I. INTRODUCTION 

Frustrated quantum spin systems on regular two- 
dimensional (2D) lattices have been the subject of intense 
interest in recent yearsi^"— They exhibit a wide variety of 
different types of ordering and phases, even at zero tem- 
perature {T — Q). Examples include various quasiclassi- 
cal antiferromagnetic (AFM) phases (e.g., with Necl or 
columnar striped ordering), phases with quantum spiral 
ordering, valence-bond crystalline phases with nematic 
ordering, and spin-liquid phases. Of particular interest 
are the (T = 0) quantum phase transitions that can occur 
as the coupling constants in the Hamiltonian are varied, 
so that the degree of frustration between bonds compet- 
ing for various types of order is changed. The result- 
ing interplay between magnetic frustration and quantum 
fluctuations has been seen to be a very effective means to 
create (and destroy) new types of order not present in the 
classical counterparts of the models. The successful syn- 
theses of ever more quasi-2D magnetic materials, and the 
experimental investigation of their properties, has also 
served to intensify their theoretical study. The very re- 
cent prospects of being able to realise spin-lattice models 
with ultracold atoms trapped in optical lattices^ is likely 
to make even more data available about the quantum 
phase transitions in the models as the exciting possibility 
opens up in such trapped-atom experiments to tune the 
strengths of the competing magnetic bonds, and hence 
to drive the system from one phase to another. 

Since quantum fluctuations tend to be largest for the 
smallest values of the spin quantum number s, for lower 
dimensionality D of the lattice, and for the smallest co- 
ordination number z of the lattice, spin-i models on 
the (hexagonal or) honeycomb lattice play a special role 
for D ^ 2, since the honeycomb lattice has the low- 
est z (= 3) of all regular 2D lattices. Thus, for ex- 
ample, one of the few exactly solvable models on the 
honeycomb lattice, namely the Kitaev model,— has been 
shown to sustain a spin-liquid phase. Clearly, the honey- 
comb lattice is also relevant to the study of graphene, 
for which much of the physics may be describable in 
terms of Hubbard-like models on this latticei^ Evidence 



has also been found from quantum Monte Carlo (QMC) 
studies^ that quantum fluctuations are sufficiently strong 
to establish an insulating spin-liquid phase between the 
nonmagnetic metallic phase and the antiferromagnetic 
(AFM) Mott insulator phase, when the Coulomb repul- 
sion parameter U becomes moderately strong. For large 
values of U the latter phase corresponds to the pure 
Heisenberg antiferromagnet (HAFM) on the bipartite 
honeycomb lattice, whose GS phase exhibits Neel LRO. 
However, higher-order terms in the t/U expansion of the 
Hubbard model may lead to frustrating exchange cou- 
plings in the corresponding spin-lattice limiting model, 
in which the HAFM with nearest-neighbor (NN) ex- 
change couplings is the leading term in the large- [/ ex- 
pansion. Frustration is easily incorporated via competing 
next-nearest-neighbor (NNN) and maybe also next-next- 
nearest-neighbor (NNNN) bonds. Recent calculations of 
the low-dimensional material /3-CU2V2O7 also show that 
its magnetic properties can be described in terms of a 
spin-i model on a distorted honeycomb latticci^ 

For all these and other reasons, frustrated spin-i 
Heisenberg models on the honeycomb lattice, including 
couplings Ji, J2, and J3 up to third nearest neighbors, 
have been extensively studied using a variety of theoret- 
ical tools.—"— The Hamiltonian for the system is 
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where i runs over all lattice sites on the lattice, and j 
runs over all NN sites, k over all NNN sites, and I over 
all NNNN sites to i, respectively, counting each bond 
once and once only. Each site i of the lattice carries a 
particle with spin operator s; and spin quantum number 
s = i . The lattice and the exchange bonds are illustrated 
in Fig. m 

The solution of the classical version of the model (i.e., 
when s — > 00) is itself richj^iiS For the AFM version of the 
model (i.e., when Ji > 0) there are six different ground- 
state (GS) phases, comprising three different collincar 
AFM phases, the ferromagnetic (FM) state, and two dif- 
ferent hehcal phases (and see, e.g., Fig. 2 of Ref.[l3)- The 
AFM phases are the Neel phase (N) shown in Fig. [1] and 
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FIG. 1: (Color online) The Neel state and the bonds of the Ji- 
J2-J3, honeycomb model; the arrows represent spins located 
on lattice sites •. 



the so-called striped (S) and anti-Ncel (aN) phases. The 
S, aN, and N states have, respectively, 1, 2, and all 3 NN 
spins to a given spin antiparallel to it. Equivalently, if 
we consider the sites of the honeycomb lattice to form a 
set of parallel sawtooth (or zigzag) chains (in any one of 
the three equivalent directions) , the S state comprises al- 
ternating FM chains, while the aN state comprises AFM 
chains in which NN spins on adjacent chains arc paral- 
lel. There are actually infinite manifolds of non-coplanar 
states degenerate in energy with each of the S and aN 
states at T = 0, but both thermal and quantum fluctua- 
tions select the collinear configurations At the classical 
level there is an exact symmetry between the GS phase 
diagrams of the AFM (Ji > 0) and FM (Ji < 0) models, 
whereby one maps into the other under the interchanges 
Ji ^ —Ji, J3 ^ —J3, and Si ^ —Si, for sites i belonging 
to one of the two equivalent triangular sublattices of the 
honeycomb lattice (c.f.. Figs. 2 and 3 of Ref . [To|) . When 
all three bonds are AFM in nature, the three possible GS 
phases (viz., the N, S, and one of the helical phases) meet 
at a tricritical point at J3 = J2 = Ji/2. 

The line J3 = J2 (= aJi) is thus of special interest, and 
we henceforth restrict ourselves to this situation where 
the NNN and NNNN bonds have equal strength for the 
remainder of the paper. There are then 4 GS classical 
phases. For the AFM case (with Ji > 0) we have: (a) 
the AFM S state when a > 5; (b) the AFM N state when 
— 1 < a < ^; and (c) the FM state when a < — 1. For 
the FM case (with Ji < 0) wc have: (a) the FM state 
when a > — j^; (b) a spiral state when — i < a < — -tj; 
and (c) the AFM S state when a < — ^ . The N state has 
classical first-order phase transitions to both the S state 
and the FM state, while the spiral state has continuous 
transitions to both the S state and the FM state. 

In two previous paper a^^'^*^ we have applied the cou- 
pled cluster method (COM) to the s = ^ version of 
this model (with J3 = J2). For the AFM case (with 
Ji > 0) we foundi^ that the direct classical first-order 
phase transition between the two (N and S) AFM states 
at a^i~^ = 0.5 is changed for the quantum s = | model 
into two separate transitions, so that a Neel-ordered (N) 
phase exists for a < ~ 0.47 and an AFM stripe- 
ordered (S) phase exists for a > ~ 0.60. In between, 
for Q!ci < a < , we found a paramagnetic GS phase 
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FIG. 2: (Color online) CCM LSUBm results for the GS en- 
ergy per spin, E/N , (Ji = 1 and J3 = J2) based on the 
Neel state (left thick curves) and striped state (right thin 
curves) as model states. We show calculated results for 
m — {6, 8, 10, 12} and the extrapolated LSUBoo result (see 



text). The exact FM result, 
shown. 



E'^^/N = 3(1 + 3J2)/8, is also 



with plaquette valence-bond crystalline (PVBC) ordering 
that has no classical counterpart. We further found that 
the quantum critical point (QCP) at appears to be 
first order, while that at is continuous. Since the N 
and PVBC phases break different symmetries we argued 
that our results favor the dcconfincment scenario^^ for 
the latter transition. For the FM case (with Ji < 0) we 
found^S that the two classical transitions from the spiral 
phase to the S phase at acf"^ = — 0.2 and from the spiral 



sp-S 

\i 

phase to the FM phase at a^^"^^ = —0.1 are changed 
for the quantum s = ^ model into one of two scenarios, 
namely, either a direct first-order transition between the 
the AFM S state and the FM state at ~ -0.11, or 
there exists an intervening phase between them in the 
very narrow range —0.12 ^ a —0.10. 

In order to complete the phase diagram of the spin-i 
model there remains to investigate the quantum analog 



of the classical first-order transition at a 



N-FM 
:1 



-1 be- 



tween the AFM N state and the FM state for the AFM 
case (with Ji > 0), and that is the purpose of this paper. 
We shall again use the CCM to do so. The two states for 
the classical model have respective energies per spin given 



J2 



Hence at the phai 
energy per spin is Ed/N 



and£;™/7V= |52(Ji+3J2). 

,N-FM 



Hence at the phase transition point a^^' = — 1 the GS 

— I if we take s = | and 
Ji = +1 to set the energy scale. In all that follows we 
take Ji = 1. 



The CCM (see, e.g., Refs. I23l - l25l and references cited 
therein) is one of the most powerful and most versatile 
modern techniques in quantum many-body theory. It 
has been applied to ma ny qua ntum magnets with huge 

success (see Refs. IHIlOiillli-lli and references cited 

therein). The interested reader can find details of the 
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FIG. 3: (Color online) COM LSUBm results for the GS 
order parameter, M, (Ji = 1 and J3 = J2) based on 
the Neel state (left thick curves) and striped state (right 
thin curves) as model states. We show calculated results 
for m = {6, 8, 10, 12} and the extrapolated LSUBoo(l) and 
LSUBcxd(2) results (see text). 



CCM in the references cited, and we do not elaborate 
here. We note only that it is a size-extensive method 
that provides results from the outset in the infinite-lattice 
limit (N — > 00). The method requires us to provide 
a model (or reference) state, with respect to which the 
quantum correlations are expressed. Here we simply use 
the N state shown in Fig.[Tl although for comparison pur- 
poses we also display below results obtained prcviouslyi^ 
based on the S state. As before, we use the well-tested 
localized lattice-animal-based subsystem (LSUBm) trun- 
cation scheme in which all multispin correlations are re- 
tained in the CCM correlation operators over all distinct 
locales on the lattice defined by m or fewer contiguous 
sites. The method of solving for higher orders of LSUBm 
approximations is discussed in detail in Ref. [2^ . 

The number of independent fundamental clusters in- 
creases rapidly with the LSUBm truncation index m. 
Hence, it is essential to employ parallel processing tech- 
niques and supcrcomputing resources for larger values 
of To obtain results in the (exact) m — > 00 limit, 
we need to extrapolate the raw LSUBm data. Since the 
hexagon is a fundamental element of the honeycomb lat- 
tice we use LSUBm data only with m > 6. For the 
GS energy per spin we employ the usual and well-tested 
scheme, E{m)/N = aq + ai?Ti~^ + a2m~'*. For the mag- 
netic order parameter (or average onsitc magnetization), 
M. different schemes have been used for different sit- 
uations. For models with no or only little frustration 
a well-tested scheme is M(m) = 5o + bim^^ + 62m~^, 
whereas a more appropriate scheme for highly frustrated 
models, especially those showing a GS quantum phase 
transition, is M{m) = cq + cim~^/^ -|- C2m~^/^. We 
henceforth refer to these latter two schemes for M as 
LSUBoo(l) and LSUBoo(2), respectively. AU of the the 
extrapolations shown below are based on LSUBm results 



FIG. 4: (Color online) The phase diagram of the spin-| Ji- 
Ji-J'i honeycomb model in the J1-J2 plane, for the case J3 = 
J2. The continuous transition between the Neel and PVBC 
phases at J2/J1 = a = ~ 0.47 is shown by a broken 
line, while the first-order transition between the PVBC and 
striped phases at ~ 0.60 is shown by a solid line. The 
transition between the striped and FM phases is either a first- 
order one at ~ —0.11, shown by a solid line, or occurs 
via an intermediate phase, probably with noncollinear spiral 
order, which exists in the region —0.12 < a < —0.10. The 
first-order transition between the FM and the Neel phases at 
Qfc4 ~ —1.17 is shown by a solid line. 



with m = {6,8,10,12}. 

II. RESULTS 

In Fig. [2] we show our CCM results for the GS energy 
per spin, E/N. They are evidently very well converged 
for all values of J2 shown. There is a clear energy cross- 
ing of the FM and (extrapolated) N energy curves at a 
value ~ -1.17, with E/N « -0.941, which is di- 
rect evidence of a first-order phase transition, just as in 
the classical case where it occurs at a^i^^'^ = — 1 with 
E^i/N = -0.75. 

We note that the individual LSUBm energy curves 
based on the Neel model state terminate at some lower 
critical value, Q;LSUBm .^^jji^h itself depends on the in- 
dex m. These termination points, below which no real 
solutions to the coupled CCM equations exist, arc them- 
selves a reflection of the actual QCP at . For example, 
Q,LSUBi2 ^ —1.38. In Fig. [3] we show the corresponding 
results for the GS magnetic order parameter, M. We 
observe that the behavior of M on the Neel side near 
the QCP at is quite smooth, with the only indica- 
tion of the transition to the FM state being the down- 
turn very near the a]f^^^"^ termination points. This is in 
sharp contrast to the behavior at the other end near the 
QCP of the N state with the PVBC state at 0;^ where 
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M — >■ 0. Clearly the best estimate for ofci comes from 
the LSUBoo(2) extrapolation, whereas the best estimate 
for M at the pure honeycomb HAFM point (i.e., when 
J3 = J2 = 0) comes from the LSUBoo(l) extrapolation, 
which gives M = 0.272 ± 0.002, in excellent agreement 
with the value M = 0.2677 ± 0.0006 from a QMC sim- 
ulation of lattices up to size N = 2048 1^ Figure |3] also 
clearly shows the corresponding transition at between 
the S state and the PVBC state, where again M — > on 
the striped phase side. 

III. SUMMARY AND CONCLUSIONS 

In this and previous papers ; ^^1^" we have studied the 
spin-i J1-J2-J3 Heisenberg model, with J2 = J3, on the 
honeycomb lattice, using the CCM. In the present paper 
we have concentrated on completing the phase diagram. 
In particular we find that the classical direct first-order 
phase transition for the AFM case (where Ji > 0) be- 
tween the AFM Neel-ordered phase and the FM phase is 
preserved for the quantum spin-i model, but now occurs 



at a QCP, ~ — 1.17± 0.01, compared to the classi- 
cal value a^"^"^ = — 1. Thus quantum fluctuations act 
to stabilize the coUinear AFM order at the expense of 
the FM order, to higher values of frustration than in the 
classical case, as has also been observed in the FM ver- 
sion of the spin-i J1-J2 model on the square latticei^ 
We find no evidence that quantum fluctuations permit 
an intervening state with no classical counterpart, unlike 
the case of the transition between the two (Neel-ordered 
and stripe-ordered) AFM states which occurs as a direct 

first-order phase transition at a!!t~^ = 0-5 for the classi- 

.1 

cal model, but occurs in the quantum spin- 2 model via 
the intermediate PVBC phase. 

Our results from this and our previous CCM 
studiesi^i2£ are summarised in the complete phase dia- 
gram for the model shown in Fig. 
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